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ABSTRACT 

Wc model the evolution of the mean galaxy occupation of dark-matter halos over the 
range O.K z < 1.3, using the data from the VIMOS-VLT Deep Survey (VVDS). The 
galaxy projected correlation function Wp{rp) was computed for a set of luminosity- 
limited subsamples and fits to its shape were obtained using two variants of Halo 
Occupation Distribution models. These provide us with a set of best-fitting parame- 
ters, from which we obtain the average mass of a halo and average number of galaxies 
per halo. We find that after accounting for the evolution in luminosity and assuming 
that we are largely following the same population, the underlying dark matter halo 
shows a growth in mass with decreasing redshift as expected in a hierarchical structure 
formation scenario. Using two different HOD models, we see that the halo mass grows 
by 90 % over the redshift interval z=[0.5,1.0]. This is the first time the evolution in 
halo mass at high redshifts has been obtained from a single data survey and it follows 
the simple form seen in N-body simulations with M{z) = Moe~^^, and /? = 1.3 ±0.30. 
This provides evidence for a rapid accretion phase of massive halos having a present- 
day mass Mq ~ 10^^-^h~^MQ, with a m > O.IMq merger event occuring between 
redshifts of 0.5 and 1.0. Futhermore, we find that more luminous galaxies are found to 
occupy more massive halos irrespectively of the redshift. Finally, the average number 
of galaxies per halo shows little increase from redshift 1.0 to z~ 0.5, with a sharp 
increase by a factor ^3 from Z'^ 0.5 to Z'^ 0.1, likely due to the dynamical friction of 
subhalos within their host halos. 

Key words: surveys - methods: statistical - galaxies: high redshift - large scale 
structure of universe 



1 INTRODUCTION 



* Based on data obtained with the European Southern Obser- 
vatory Very Large Telescope, Paranal, Chile, program 070. A- 
9007(A), and on data obtained at the Canada- France- Hawaii- 
Telescope, operated by the CNRS of France, CNRC in Canada 
and the University of Hawaii, 
f E-mail: abbas@oato.inaf.it 



The correlation function of galaxies is a simple yet powerful 
tool that allows one to constrain cosmological parameters 
and models of galaxy formation. Furthermore, with the help 
of high redshift surveys, the evolution in the clustering of 
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galaxies allows for a better discrimination b etween theoret i- 
cal models degenerate at the present epoch l|Peacocklll997l ). 

Until recently, the galaxy corre lation function had been 
thought to follow a power law (|Totsuii fc Kiharal 1 19691 : 
|Peeblej[T97i : iGott fc Turneilll979l '). Subsequently, a depar- 
ture from the power law on small scales (of the order 1 to 
a few Mpc/h) in the galaxy correlation func tion was no- 
ticed in pioneering surveys of the eighties, e.g. (|Guzzo et al.l 
and has now been fully confirmed by many large 
and deep galaxy surveys. These consist of various surveys 
at low and intermediate rcdshi fts (z ^ 1.5) s uch as the 
SPSS ([Cp nnolly ct al. 2002; Z ehavi et all |2004|). 2dFGRS 
jMagliocchetti fc Porcianil |20o3), VV DS jLe Fevre et all 



l2005bl : iPollo et al.ll2006lC cOMBO-17 (|Phleps et al.ll2006h . 
DEEP2(| Coil et al.1 bOOd ) and for lyman break galaxies 



(LBGs) at hi gh redshifts in the SXDS and GOODS surveys 
jOuchi et al.i r2005: Lee ct al. 200ii). 

Earliest measures showing a non-power law for the 
galaxy correlation function were difficult to interpret as they 
relied heavily up on the angular correl ation function, u]{9), 
at low redshifts (|Connollv et al.ll2002l ). This meant an in- 
tegration over a wide range of galaxy luminosities and red- 
shifts as well as complicated co rrelations between the statis- 
tical errors l|Zehavi et all 120041 '). This was followed by mea- 
surements of the projected correlation function, lOp(rp), and 
uj{6) in larger and deeper surveys. It has been seen that 
the deviation from a power law becomes more pronounced 
for bright galajcy samples with L > Lt (ICoil et al.l l2006j: 



Polio et allbooel ) and for LBGs at high redshifts (|Lee et al.1 
20061 ). Likewise, in SPH simulations the luminous and more 



strongly clustered galaxies show a similar behaviour where 
the 'kink' is clearer tha n in the case of the full galaxy sample 
l|Weinberg et al.ll2004 ). 

Interestingly enough, the correlation function for dark 
matter in N-body simulations is well known to be non- 
adherent to a power law (Jenkins et al. 1998, Kauffmann 
et al. 1999, Cooray & Sheth 2002 and references therein). 
The natural question that arises is how biased is the galaxy 
distribution with respect to the underlying matter distri- 
bution? The overall shape of the dark matter correlation 
function is mostly unaffected as one goes to higher redshifts 
as seen in SPH si mulations (^Weinb erg et al.l |2004| ) and N- 
body simulations (jjenkins e t al. 19 9^. This is different to 
what is seen for hig h-z galaxies (|Ouchi et al.ll2005l : lLee et al.l 
l2006l : IZheng et al.ll2007i where the so-called 'break' is more 
prominent implying that the biasing of the galaxy distribu- 
tion on large scales, spatial exclusion of dark matter halos 
on small scales, alongwith a host of other complex physical 
processes, such as dynamical friction, feedback from super- 
novae, ram-pressure stripping etc. conspire in a non-trivial 
way to produce differences in the galaxy correlation function 
at different redshifts. 

The break in the power law can be physically inter- 
preted in the language of the halo model (see Cooray & 
Sheth 2002 for a detailed review) as the transition be- 
tween two scales - small scales lying within the halo to 
those larger than the halo. It is only natural to use a halo- 
based prescription where g alaxies form by the co oling of gas 
within dark matter halos jWhite fc Reed 1 19781 ). which are 
bound, virialized clumps of dark matter that a re roughly 200 
times the background density at that time (|Gunn fc GottI 
1 19721 ). The galaxies occupy dark matter halos following a 



HOD (halo occupation distribution) model. In turn the 
HOD fully describes the bias in the distribution of galax- 
ies with respect to the und erlying dark matter distribution 
((Berlind fc Weinbergj[20o3 ). 

The motivation for HOD based models axose when 
it was noticed that the clustering of galaxies could be 
reproduced by populating halos in semi-analytic models 
with galaxies following a particular probability distribu- 
tion citekau99,ben00,ber02. Furthermore, it has been seen 
that without the help of a proper halo based description 
the strong clustering of red galaxies (at z ~ 3) can be ex- 
plained by high and unrealistic (anywhere in the range of 
70-200 galaxies per halo) occupation numbers to match the 
observed number density a nd strong cl ustering of a small 
number of high mass halos (|Zheng||2004l ). 

Recently, several groups have studied the galaxy corre- 
lation i n light of the HOD mo dels (Magliocchctti & Porcianil 
20031: Ivan den Bosch et al.1 [20031: iHamana et al.1 |200l 



Zehavi et al 



Conrov et al.l 



2005 



Ouchi et al 



20061 : IZheng et al 



20051 : iPhleps et al.1 l2006l : 
20071 ). Most of these works 



have mainly concentrated on obtaining best-fit HOD pa- 
rameters and consequently the evolution in the HOD and 
information on the underlying dark matter distribution. In 
some cases, data from different surveys havin g different se- 
lections were used to s tudy the evolution (|Conrov et al.l 
I2OO6I : IZheng et ahi |2007') The work in this paper comple- 
ments and extends these analyses by studying the evolution 
in the HOD over a redshift range z ~ 0.1 - 1.3 for a variety 
of luminosity-limited samples but always for data from the 
same survey, the VIMOS VLT Deep Survey (VVDS). We 
also study two different HOD models in order to obtain a 
better understanding on the degeneracies between the vari- 
ous parameters. 

The paper is organized as follows. In section[2]we briefly 
describe the data used from the VVDS survey. This section 
is followed by Section |3] giving an outline of the theoret- 
ical framework. Section |4] describes the fltting procedure 
alongwith the best-fit parameters obtained for the differ- 
ent models. The estimates for the average halo mass and 
number of galaxies per halo (galaxy weighted) for the differ- 
ent luminosity-threshold samples are also presented. Finally, 
Section [S] wraps up with a discussion and conclusion of the 
results. 

Throughout we will assume a flat ACDM model for 
which {flo,h,as) = (0.3,0.7,0.9) at z = 0. Here Qq is the 
density in units of critical density today, h is the Hubble 
constant today in units of 100 km s~^ Mpc~^, and ag de- 
scribes the rms fluctuations of the initial field evolved to the 
present time using linear theory and smoothed with a top- 
hat filter of radius 8 Mpc/h. All absolute magnitudes are in 
the AB system. 



2 THE VVDS DATA AND ITS ANALYSIS 
2.1 Description of the data 

The data used in this analysis comes from the First Epoch 
VVDS (VIMOS- VLT Deep Survey) 0226-04 "Deep" field 
(hereafter VVDS-Deep). A complete description of the data, 
survey strategy, data reduction and primary goals can be 
found in Le Fevre et al. (2005a). Here we will simply give a 
short description of the data used. 
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Table 1. Properties of the different subsamples 
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Figure 1. Mb for the VVDS galaxies at different redshifts. The 
solid colored lines denote the various subsamples selected with 
properties mentioned in Table [T] The dashed line shows the ob- 



The VVDS-Deep sample is magnitude limited in the 
Iab band with 17.5 ^ Iab ^ 24 and covers an area of 

0. 49 deg^ without any color or shape restrictions imposed. 
The spectroscopic observations were taken with the Visi- 
ble Multiple-Object Spectrograph (VIMOS, Le Fevre et al. 
2003) at the ESO VLT, whereas the Virm os Deep Imag- 
ing s urvey (VDIS) BVRI photometric data l|Le Fevre et al.l 

was obtained with the wide field 12k mosaic cam- 
era at the CFHT (Canada-France-Hawaii telescope) and is 
complete and free from surface brightness selection effects 
l|McCracken et al.ll2003l ). 

The sample contains 6582 galaxies with secure redshifts, 

1. e. known at a confidence level ^ 80%. These galaxies have 



a mean redshift of z ~ 0.83. Fig. \T\ shows the absolute 
magnitude of these galaxies as a function of redshift along- 
with the different luminosity-threshold subsamples selected, 
where the luminosity threshold is assumed to evolve accord- 
ing to the relation Afs (z) = —1.15z + Mb{z — 0). The factor 
of '- 1.15' arises from the redshift evolution of the character- 
istic absolute magnitude, Mg , of galaxies as measured in the 
luminosity function. This value has been determined using 
the luminosity function measurements obtained within the 
same sample by Ilbert et al. (2005). 

An evolving luminosity threshold needs to be taken into 
consideration when comparing samples at different epochs, 
as it provides us with statistically similar samples at differ- 
ent redshifts, having similar evolved luminosities. Assuming 
that the global evolution of galaxies has as a main conse- 
quence to increase the global luminosity of galaxies, we fol- 
low the evolution of galaxies with similar properties on av- 
erage. This falls within the boundaries of standard practice 
of galaxy evolution studies. 

As our subsamples are nearly volume complete, and as 
we are using all types of galaxies together, we may follow 
the global increase in the halo mass of an average galaxy. 
However, we do recognize that this way of selecting galaxies 
does not garantee to follow the exact same population with 
cosmic time. Unfortunately, there is no single prescription 
enabling to tag galaxies and exactly follow their precursors 
/ descendants. Indeed, if this were possible it would be the 
solution to galaxy evolution. To try to quantify the impact 
of our selection on the average halo mass, we have used the 
Millennium simulation. This will be further discussed in the 
next section. 

Samples using a similar type of selection, i.e. using lu- 
minosity thresholds, ha ve been extensively studied within a 
theoretical framework (Zehavi et al. 20051 : ICoil et aL I I2OO6I : 



IConrov et al] l2006l : IZheng et all 120071 ). The corresponding 



HOD parametrisation requires fewer parameters to be fit- 
ted as compared to differential, luminosity binned samples. 
However, this means that one is biased towards increasingly 
brighter galaxies at higher redshifts, simply due to the fact 
that the sample is selected in apparent magnitude. As can be 
seen in Fig. [TJ there is a change of about 5 in the B-band ab- 
solute magnitude over the redshift bin z ~ [0.1, 1.5]. Table [1] 
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shows the various properties of the subsamples alongwith 
their number densities. The galaxy number densities were 
computed by integrating the luminosity functions, derived 
by Ilbert et al. 2005 on the same galaxy sample and parame- 
terized using Schechter functions. The evolution in the best- 
fit Schechter parameters, M* , (j>* and a was taken into con- 
sideration, thereby accounting for luminosity evolution. We 
estimated the errors on the number densities by propagat- 
ing the errors on the Schechter parameters. M^^'^"'^ denotes 
the evolving absolute magnitude threshold at the highest 
redshift. For each M^g^'^"^ two samples were obtained, one 
at low redshift and another at higher redshift with brighter 
galaxies selected due to the evolving selection cut. The two 
samples overlap slightly in redshift in order to maximise the 
number of objects. Mg^ and M^^^ are the absolute mag- 
nitudes of the evolving cut at the lower and higher redshift 
limits respectively of the redshift range. 

2.2 Studies with simulations 

Ideally, we would like to follow statistically the same galaxy 
population with time in order to study the growth of the un- 
derlying dark matter halo mass. This is a tricky issue as it is 
difficult if not impossible to know the exact progenitors of a 
descendant galaxy population and how to select them. Ifow- 
ever, by taking care to follow the exact population mix down 
to a fixed absolute luminosity, we can minimize the bias in 
the average halo mass due to the presence/lack of faint or 
bright galaxies. As mentioned above, this is made possible 
by accurate measurements of galaxy evolution. In order to 
tackle this issue we use the Milli-Millennium simulation of 
galax ies having 270^ part icles in a box of 62.5Mpc/h on its 
side (jSpringel et al.ll2005l ). The simulations retain informa- 
tion on the progenitor trees of galaxies making it possible 
to make a comparison with volume-limited samples. 

Let us select a galax;y sample at high redshift chosen 
with a luminosity cut-off in the B-band similar to what is 
done in the VVDS data, hereafter called the 'parent sam- 
ple'. This sample is then evolved into two samples at lower 
redshift, a 'simulated sample' at the lower redshift, having 
a luminosity cut-off that is evolved and fainter (again sim- 
ilar to what we did in the data), and another sample that 
contains all the descendant galaxies at the same lower red- 
shift (hereafter the 'descendant sample'). Doing a galaxy to 
galax;y match between the two samples would tell us how 
many galaxies in the simulated data set are actual descen- 
dants and therefore the same population followed through 
time and the effects on the underlying average halo mass. Ta- 
ble [2] shows the Millennium samples selected having roughly 
the same mean redshift and mean absolute magnitude, AIb, 
as in the VVDS data sample of Table [1] 

We find that at worst 77% of the M < -20 and at best 
87% of the A/ < —18 simulated luminosity threshold sample 
are actual descendants at lower redshifts. From Table [2] we 
can also study the effect of the selections on the underlying 
average halo mass. It can be seen that typically the under- 
lying halos in the descendant sample are heavier than those 
in the simulated sample. After taking a closer look at the 
descendant sample we noted that even though there are a 
larger number of fainter galax;ies, there are also slightly more 
bright galaxies that lead to a slightly higher average magni- 
tude. The combination of faint satellite galaxies residing in 



massive halos and fewer galaxies of intermediate luminosity, 
likely lead to a descendant sample with more massive halos 
and slightly brighter galaxies on average than the simulated 
sample. This possibly causes a lower overlap between the 
simulated and descendant samples for the brightest samples. 

The Millennium simulation shows that a growth in halo 
mass detected in the data would be underestimated with re- 
spect to what could be seen ideally. The underestimation in 
mass is of the order of roughly 10 %, and therefore a measure 
in the growth of mass of a halo can be mainly attributed to 
the hierarchical formation of structure and not due to the 
typology of the selection (taking into consideration the high 
overlap between the simulated and descendant samples). 



2.3 The correlation function 

The redshift-space correlation functions for the different lu- 
minosity threshold samples have been computed via the 
Landy & Szalay (1993) estimator: 

. Nr{Nr - 1) GG{rp,n) Ne - 1 GR{rp,7T) 
^^''P'''' NciNo^l) RR{rp,7v) No RR{r^,n)^ 

(1) 

where No and A''^; are respectively the total number of 
galaxies and randomly distributed points in the same sur- 
vey. GG{rp, n) is the number of distinct galaxy-galaxy pairs 
with separations lying in the interval {tt,tt + dn) in the ra- 
dial direction and {rp,rp + dvp) perpendicular to the line of 
sight. Likewise, RR{rp,-K) and GR{rp,n) are the number of 
random-random pairs and galaxy-random pairs respectively 
in the same interval. 

In order to avoid redshift space distortions, ^{rp, n) has 
been integrated along the line of sight to obtain the pro- 
jected correlation function ([Davis fc Peebles! [l983l l: 

i^p{rp) = 2 / ^(r-p,7r)d7r 

■^0^ (2) 

= 2^"c[\AiT7]dy, 

where ^(r) is the real space correlation function with r = 
+ y^. The measu rements using the same sample impose 
a sim ilar upper limit (jPollo et al ] |2005l . I2OO61 : .Meneux et all 
I2OO6I I. Polio et al. (2005) found that LOp{rp) is quite insensi- 
tive to TVmax in the range of 15Mpc/h < nmax < 25Mpc/h 
for rp < lOMpc/h. Too small a value for this limit would 
cause an underestimation of the small-scale power, and too 
large a value would introduce noise. After several exper- 
iments, the optimal value of ■Kmax='20 Mpc/h has been 
adopted. 

The errors have been estimated using bootstrap resam- 
pling of the data, which consists of computing the variance of 
Wp{rp) in Nreai bootstrap realizations of the sample. Each 
realization is obtained by randomly selecting a subset of 
galaxies from the data sample allowing for repetitions. A 
correction factor is then applied to account for the under- 
estimation of the errors obtained using this technique. This 
correction factor has been calibrated on mock samples to 
match the ensemble error (accounting for cosmic variance) 
of the simulated mock samples (see Polio et al. 2006). 
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Table 2. Simulated data sets vs. descendants 
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Figure 2. The correlation function for the various luminosity threshold samples. The symbols and error bars denote the measurements 
from the VVDS-Deep. The lines present the best-fit halo model for Z model of Zehavi et al. (2005) in the left panel and for TWZZ model 
of Tinker et al. (2005) in the right panel. 



3 ANALYTICAL MODELING 
3.1 The Halo model 

The analytical model is based on the halo model (see Cooray 
& Sheth 2002, for a review), here we will briefly mention the 
main ingredients. All mass is assumed to be bound up into 
dark matter halos having a range of masses which in turn 
host galaxies. 

In this model the power spectrum, P(k), and/or cor- 
relation function, ^(r), of the galaxies (which are fourier 



transform pairs) can be written as the sum of two terms. 
One that dominates on non-linear scales -smaller than the 
size of a halo, and the other term becoming significant on 
larger linear scales, known as the 1-halo and 2-halo terms 
respectively. The 1 halo term arises from pairs of galaxies 
lying within the same halo, whereas pairs of galaxies lying 
in different halos contribute to the 2 halo term. In fourier 
space this can be written as. 
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P{k) = P^n{k) + P2h{k), 



(3) 



where, 



— — / dmn(m,, z){Ng{Ng — \)\m) \u{k\rn)\ , 
^gai Jo 

- — / dmn{m, z){Ng\m) b{m, z) \u{k\m 

ngal Jq 

X Plin{ky z') , 

(4) 



and Ugai is the number density of galaxies: 



I'gal — 



dmn{m, z){Ng\m). 



(5) 



Here, {Ng\m) is the average number of galaxies occupying a 
halo of mass m, u{k\m) is the halo density profile in fourier 
space, n(m, z) is the number density of halos of mass m, 
b{m, z) is the bias factor which describes the strength of 
halo clustering, and Pun (fc, z) is the power spectrum of the 
mass in linear theory all at a given redshift z. The upper 
limit of integration, Mmax, approximately accounts for the 
halo exclusion efi'ect (different halos cannot overlap) by sup- 
pressing the 2-halo term at small scales. Following Zehavi et 
al. (2004), Mmax is the mass of the halo with virial radius 
r/2. One can also calculate the one-halo term for the correla- 
tion function exactly in real space, which is the approach we 
have taken. For more details we refer the reader to Berlind 
& Weinberg (2002). The two-halo term has been computed 
in k-space and then fourier transformed to obtain the cor- 
relation function. Then the projected correlation function is 
obtained as in Equation [2] Similarly to the data the upper 
limit is chosen to be finite {rmax='20 Mpc//i) in order to 
avoid noise caused by uncorrelated distant pairs. 

We assume that the density profiles of halos have 
the form described by Navarro et al. (1997), with a halo 
concentration parameter c{M) = 11(M/M«)~'^ '^^ to ac- 
count for the definition of h alos as spheres en closing 200 
times the background density jZehavi et al.|[2005t ) and where 
o"(M, , z) — Ssc (o" and Ssc are defined below). The halo abun- 
dances and clustering are described by the Sheth & Tormen 
(1999) parameterization: 

771 

— n{m, 2)d77i = f{i')di' 

= -7r\ -^""^P — ^ + ) ] 



6(777, z) = 1 + 



av^ - 1 



2 

2p/5, 



1 + {au'^y 



a(m, z) ' 



allowing us to write the background density as, 



P= d777 77(777, 2) 777. 



(6) 



(7) 



with Ssc being the critical density required for spherical col- 
lapse, extrapolated using linear theory to the present time 
(= 1.686, ignoring the weak cosmological dependence), a ~ 
0.71, p = 0.3, and A ~ 0.322. a{m, z) is the rms value of the 



initial fluctuation field when smoothed with a top hat filter 
and extrapolated to the present time using linear theory. 



(7(777, z) — aim, 0) 



G{z) 
G(0) 



(8) 



where G is the growth factor from Carroll, Press & Turner 
2 (1992). 



3.2 The HOD models 

We will consider two similar HOD models. The first one is 
based on the model used in Zehavi et al. (2005) (hereafter Z 
model) to compare to the SDSS data, and is motivated by 
Kravtsov et al. (2004). 



(iV,|777> = l + 







Ml 



for m > Mmin 



otherwise 



(9) 



The second model was proposed by Tinker et al. (2005) 
(hereafter TWZZ model) and is given by, 



{Ng\m) = 1 + ;^exp 
= 



Mc 



for m > M„ 



otherwise 



(10) 



where Mmin is the minimum mass for a halo to host one cen- 
tral galaxy, and Mi is the mass of a halo hosting on average 
one satellite galaxy. The T' represents one central galaxy 
placed at the center of mass of the parent halo, and the satel- 
lite galaxies follow the underlying dark matter distribution. 
TWZZ model has been used to study the HOD for a range 
of redshifts and number densities and found to give results 
for the correlation function in good agr eement with variou s 
redshift surveys in the range z = [0-5] (|Conrov et al ] |2006l ). 

Our purpose here will be to obtain the best-fit HOD 
parameters for the two models, and compare the number- 
weighted halo masses and number of galaxies. We have de- 
cided to use these models in order to keep things simple and 
easy to interpret. Based on the statistics of the sample and 
the number of data points used, it is best to use HOD mod- 
els with a minimal number of free parameters adapted to 
the science case at hand, i.e. for this paper, the average halo 
masses and number galaxies mentioned below. The main re- 
sults that are obtained should remain essentially the same 
irrespective of the HOD model chosen. 

A complementary approach would be to use modeling 
based on conditional luminosity functions (CLF) or condi- 
tional occupation numbers (CON) (e.g. van den Bosch et al. 
2003, Yang et al. 2003, Cooray 2006). The present attempt 
is a first at using a large number of galaxy spectra at high 
redshift to study the evolutionary behaviour of a few prop- 
erties pertaining to the galaxy and dark matter distribution. 
These 'few' properties are certainly not comprehensive, and 
this work can be seen as a starting point for more studies 
using the sample. Moreover, larger data samples from on- 
going and upcoming redshift surveys will certainly provide 
grounds for extensive studies based on CLF/CON. 
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Table 3. Results for Z model 



M 


■high 


logioMmin 


logioMi 


a 


logio < M > 


< Af > 


xVdof 


rifu 


< 


-17.77 


11.05 ±0.87 


12.73 ± 1.00 


0.99 ±0.10 


13.00 ±0.13 


1.21 ±0.01 


1.03 


37.28 


< 


-18.00 


11.21 ±0.18 


13.03 ±0.11 


1.27 ±0.08 


12.97 ±0.08 


1.12 ±0.01 


1.35 


26.29 


< 


-18.27 


11.37 ±0.16 


12.93 ±0.06 


1.11 ±0.12 


13.14 ±0.14 


1.22 ±0.01 


0.88 


20.65 


< 


-18.50 


11.63 ±0.19 


13.31 ±0.18 


1.31 ±0.09 


12.98 ±0.06 


1.11 ±0.02 


0.52 


11.16 


< 


-18.71 


11.88 ±0.08 


13.57 ±0.10 


1.33 ±0.16 


13.10 ±0.08 


1.09 ±0.03 


1.16 


6.59 


< 


-19.00 


11.76 ±0.08 


13.39 ±0.08 


1.19 ±0.24 


12.88 ±0.09 


1.11 ±0.03 


0.84 


8.53 


< 


-19.16 


12.06 ±0.12 


13.72 ±0.10 


1.43 ±0.27 


13.11 ±0.08 


1.08 ±0.05 


1.51 


4.39 


< 


-19.50 


11.92 ±0.16 


13.43 ±0.22 


1.22 ±0.25 


12.94 ±0.07 


1.12 ±0.04 


1.24 


5.99 


< 


-19.48 


12.06 ±0.07 


13.78 ±0.04 


1.66 ±0.18 


13.13 ±0.06 


1.06 ±0.05 


1.08 


4.36 


< 


-20.00 


12.35 ±0.11 


13.80 ±0.16 


1.42 ±0.41 


13.09 ±0.04 


1.08 ±0.10 


0.78 


2.17 



4 RESULTS 



4.1 Results from VVDS 



The different parameters were allowed to vary within the 
following ranges: 10 log^M^^^) < 15, 10 log{Mi) < 15, 
10 ^ log{Mcut) ^ 15, 0.5 ^ a ^ 2.0. These limits represent 
reasonable constraints on the typical mass of a dark mat- 
ter halo and the power law slope. The minimum mass for a 
halo to host one central galaxy is usually ^ 10^^/i~^M© for 
low redshi ft galaxies (Zoh avi et al. 2005.) , and LBGs at high 
redshifts l|Hamana et al.l 120041 ) . At the high mass end, the 
brightest SDSS galaxy samples have Mmm < 10^'*/i~^Mq. 
Taking into account our sampling having brightest samples 
of galaxies at high redshifts, in the hierarchical structure 
formation scenario this then represents an upper limit on 
the mass. O n the other hand, Mi has been found to be 
~ 23Mmi„ l|Zehavi et al. I l2005l ). Furthermore, power law 
slopes ^ 1.5 are considered "artificially high", which gen- 
erally dominate bri ghter samples that have fewer satellite 
galaxies on average (jConrov et al.l [20061 ). The number den- 
sity obtained using Eq. [5] was restricted to lie within 3a 
from the observed number density given in Table [1] The 
correlation functions for the different luminosity thresh- 
old samples is shown in Fig[2] alongwith the best fits for 
the two HOD mo dels obtained with the MPFIT algorithm 
l|Markwardtl l2009D that uses the Levenberg-Marquardt tech- 
nique l|More et al.lll978l ) to solve the non-linear least-squares 
problem using the full covariance matrix. 

Tables [3] & |4] show the best-fit parameters for the two 
different models obtained by a minimum chi-square esti- 
mate, the value of the reduced x^, alongwith the average 
number-weighted halo masses and number of galaxies per 
halo defined as: 



(M) 
(iV) 



/J" n(m) {Ng \m)m dm 
/o°° w(m)dm 
n{m){Ng\m)Am 
ir w(m)dm 



(11) 



The generalized chi-square estimate is obtained the 
usual way adopting: 



"■bin "bin 

2 \ ^ \ obs / \ model/ Ir ohs / \ model f 

(12) 



where nbin is the number of bins and dj is the covariance 
of the values of Wp between the ith and jth bins. 

The results obtained from both models and given in 
the tables are found to be in agreement (at least for com- 
parable parameters), which is to be expected as the models 
are similar. In the case of TWZZ model the power law ex- 
ponent is kept constant and the number of satellites has a 
smooth, exponential cut-off. For example, in both cases the 
value for the minimum mass (Mmin) is very similar if not 
the same. The power law exponent for the Z model several 
times shows values of a that are quite high (> 1.4). Arti- 
ficially high values have bee n noticed in fits to simulations 
as well (|Conrov et al.l [20061 ) and occur for galaxy samples 
at high redshifts. The 1-a error bars on the parameters are 
obtained with the MPFIT algorithm. 

Figure|3]shows the number-weighted average halo mass, 
< M >, versus the redshift. The symbols with error bars 
represent the various sub-samples selected from the VVDS. 
The error bars are obtained based on error propagation for- 
mulas. The point at the lowest redshift (z ~ 0.1) is obtained 
from the SDSS using the best-fit HOD parameters from Ze- 
havi et al. (2005). The mass in this case was calculated using 
the Z model for the luminosity threshold sample having the 
same M — M* difference as the samples at higher redshift 
in the VVDS (where the difference in the r-band for the 
SDSS has been converted to the B-band, Ilbert et al. 2005). 
It can be seen that the halo mass evolves and increases as 
one goes to lower redshifts. This is an indication of the halo 
mass growth due to the hierarchical aggregation of matter. 
We find that on average < M > increases by 90% from 
redshift ~ 1 to ~ 0.5, showing that massive halos have a 
rapid accretion phase quite late on, similar to what is ex- 
pected from N-body simulations (Wechsler et al. 2002). As 
shown in Wechsler et al. (2002) the mass growth can be 
easily characterized by the form M{z) = Moe~'^^ . The in- 
teresting comparison with the addition of low redshift SDSS 
points gives a linear minimum fit of /3 ~ 1.94 ± 0.10 for 
the Z Model and /3 ~ 2.09 ±0.04 for the TWZZ model. This 
is to be compared to the predictions of the mass accretion 
history of halos in N-b ody simulations and halos generated 
through PINOCCH IO ([Monaco et al.ll2002l : IWechsler et al.[ 
[20021 : iLi et"al]|2007l ). where (3 ~ 0.62. One can argue that 
the direct comparison between data obtained from different 
rest-frame bands can be tricky and could in part lead to a 
)] slight boost in /3, even though necessary care has been taken 
in converting to a common rest-frame band. The latter is 
reflected in the value obtained for /3, using only the VVDS 
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Figure 3. The evolution in tlie number- weighted average halo mass given by Eg llll for various luminosity threshold samples is shown 
for the Z model (left panel) and TWZZ model (right panel). The symbols with error bars are obtained from the VVDS, whereas the low 
redshift symbols at 2 = 0.1 with small error bars are from the SDSS. The mass growth can be characterized by M(z) = MQe~^^, with 
P = 1.07 ± 0.57 (for the VVDS points), /3 = 1.94 ± 0.10 (VVDS + SDSS points) in the case of the Z model, and /3 = 1.54 ± 0.13 (for the 
VVDS points), /3 = 2.09 ± 0.04 (VVDS + SDSS points) for the TWZZ model. The values of /3 =0.62 and 2.00 respectively represent the 
prediction from N-body simulations and the VVDS + SDSS samples. 




Figure 4. The evolution in the satellite fraction for various luminosity threshold samples is shown for the Z model (left panel) and 
TWZZ model (right panel). The symbols with error bars are obtained from the VVDS, whereas the low redshift symbols at z = 0.1 are 
from the SDSS. Simply for illustrative purposes, the curves show the reciprocal standard power law behaviour, y = l/(ax^), of the data, 
with a = 13.05 and b = 0.61 for the Z model, and a = 20.41 and b = 0.79 for the TWZZ model. 
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Table 4. Results for TWZZ model 



M 


high 


logioMmin 


logioMi 


logioMcut 


log{< M >) 


<N> 




nfit 


< 


-17.77 


11.00 ±0.15 


12.54 ±0.39 


13.01 ± 2.07 


13.09 ±0.03 


1.11 ± 0.01 


1.19 


38.00 


< 


-18.00 


11.45 ±0.19 


11.87 ±0.71 


14.93 ±0.62 


12.82 ±0.05 


1.01 ± 0.02 


1.33 


14.51 


< 


-18.27 


11.67 ±0.13 


13.03 ±0.22 


13.32 ±0.75 


13.12 ±0.01 


1.09 ± 0.02 


0.75 


10.06 


< 


-18.50 


11.62 ±0.08 


12.93 ±0.08 


13.26 ±0.24 


12.98 ±0.01 


1.08 ± 0.02 


0.51 


11.09 


< 


-18.71 


11.88 ±0.09 


13.22 ±0.29 


13.45 ± 0.74 


13.10 ±0.01 


1.07 ±0.04 


1.16 


6.40 


< 


-19.00 


11.74±0.11 


13.22 ±0.15 


12.62 ±0.89 


12.88 ±0.02 


1.09 ± 0.03 


0.81 


8.71 


< 


-19.16 


12.14±0.11 


13.47 ±0.19 


13.49 ±0.61 


13.12 ±0.01 


1.06 ± 0.06 


1.42 


3.67 


< 


-19.50 


11.90 ±0.16 


13.23 ±0.35 


12.73 ±0.82 


12.94 ±0.02 


1.11 ± 0.04 


1.17 


6.18 


< 


-19.48 


12.07 ±0.05 


13.15 ±0.16 


13.97 ±0.25 


13.13 ±0.01 


1.04 ±0.05 


0.97 


4.18 


< 


-20.00 


12.32 ±0.12 


13.54 ±0.45 


13.21 ±0.76 


13.09 ±0.02 


1.07 ±0.10 


0.73 


2.27 



points leads to a smaller value (~ 1.54(1.07) ± 0.13(0.57) 
for TWZZ (Z) model) albeit with larger error than that 
obtained from the extrapolation to smaller redshifts, but 
slightly more consistent with the results from simulations. 

For samples at similar redshifts we can see that the 
number-weighted average halo mass increases with the lu- 
minosity threshold of the sample reinforcing the notion that 
luminous galaxies occupy massive halos. Thi s is in agree- 
ment with results ob tained from simulations (IConrov et al.] 
l2006l ) and for LBGs (|Ouchi et al.ll2005l : iLee et al.ll20m 



Figure S] presents the evolution of the galaxy satellite 
fraction, or the average number of satellite galaxies. The 
illustrative reciprocal power law behaviour of the dataset 
shows relatively little change in the satellite fraction (always 
close to ~ 0.1 within 1 a) over the redshift range of z=[0.5- 
1.0]. Over z=[0.1-0.5] there is a sharper increase by a factor 
~3 to the local SDSS value of ~0.3. The evolution is mainly 
accentuated by the SDSS points, although the two lowest 
redshift VVDS points for the case of the Z model do hint 
towards an increase with lower redshifts. It is possible that 
the sharper upturn is once again caused by the complicated 
comparison between two different data surveys. However, 
here again care has been taken to convert to the appro- 
priate rest-frame band when making these comparisons and 
should not affect the overall trend. The increase in the satel- 
lite fraction as one goes to lower redshifts can be explained 
by the dynamical fricti on of subhalos within their host ha- 
los jConrov et al. ||2006'). Subhalos are more likely to remain 
intact within massive halos, whereas in less massive halos 
they are subject to more dynamical friction and can easily 
be destroyed. The dynamical friction becomes more/less ef- 
ficient as a function of the relative masses of subhalos to 
distinct halos. This is to be compared to recent results ob- 
tained by Zheng et al. (2007) who find that the evolution of 
the sateUite fraction follows a trend similar to what is seen 
here. 

The following Fig. [5] shows the evolution in the halo 
occupation, Ng(m), for the extreme luminosity threshold 
samples obtained from the best fit parameters for the two 
models. Evidently, the minimum mass, Mmin increases with 
the lu minosity of the sa mple as is found locally in the 
SDSS IZehavi et all 120051 ') . again demonstrating that lumi- 
nous galaxies occupy more massive halos. 



4.2 Comparison to SDSS 

In this section we will compare to results for the same HOD 
model (Z model) as used in Zehavi et al. (2005). Figure [S] 
shows the comparison between the masses of halos that have 
at least one central galaxy (Mmin) and one satellite galaxy 
(Ml) on average as a function of Lthresh/ L, , where the ratio 
is in the B-band, and Lthreah is the luminosity threshold 
given in Table [1] (Lt^resii and L* are at similar redshifts). 
Here we will try to compare results obtained at different 
redshifts. For 40% of the VVDS samples, Mmin is similar to 
the local SDSS results within the error bars, with the rest of 
the VVDS samples having higher values of Mmin- Generally, 
the VVDS samples exhibit more massive halo masses. Mi, 
required to host satellite galaxies than what is seen locally. 
The value for the power law slope, a, is mostly similar to 
that for local galaxies, with the bright intermediate redshift 
galaxies showing a higher slope. We can see that generally 
the samples with higher values for a, also have higher values 
of Ml and Mmin than present-day galaxies. 

It is interesting to note the Mi /Mmin ratio, which is 
on average ~ 45, rather high as compared to the value of 
~ 23 for the SDSS galaxies. A direct comparison and inter- 
pretation of these results is complicated as one is looking 
at two different surveys taken in different restframe bands. 
However, we can speculate that the high value of the ratio 
implies that the halo with one central galaxy needs to ac- 
crete roughly 45 times its mass in order to host a satellite 
galaxy. In other words, a halo of a given mass is likely to 
have fewer satellite galaxies at higher redshifts as opposed 
to a halo of the same mass observed locally. 



5 DISCUSSION & CONCLUSIONS 

The comparison of analytical models and data provides 
useful information of how the distribution of galaxies de- 
pends on the underlying dark matter. Subsequently, the 
best-fitting parameters obtained as a result of this compari- 
son provide physical information regarding the dark matter 
halos and galaxies. 

The size of the VVDS dataset allows one to study, with 
a unique sample, the global change in the underlying halo 
properties of an average galaxy down to z ~ 1. We attempt 
to follow the evolution in some properties of a magnitude se- 
lected sample, evolving the magnitude cut-off based on accu- 
rate measurements of galaxy evolution. We have presented 
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Figure 5. The evolution of the halo occupation for the extreme luminosity threshold samples having Mb < -18, -20. The solid line 
corresponds to the lower redshift sample in each case, and the dashed line to the brighter sample at higher redshift both having the same 
evolved luminosity threshold (see Table[l}. The solid lines correspond to 2 = 0.44, 0.67 and the dashed lines to 2 = 0.66, 1.05 respectively 
for the samples with Mg < -18, -20. The left panel depicts Z model and the right panel T^TiZ model. 



results of the fitting of analytical lialo models, incorporating 
simple HOD models with minimal number of free parame- 
ters, to data (in this case the projected 2-point correlation 
function) from the VVDS survey. This allowed us to study 
the evolution of the average number weighted halo mass and 
satellite fraction. 

On different scales there are contributions from central- 
satellite, satellite-satellite, and central-central pairs of galax- 
ies to the correlation function thereby providing constraints 
on the evolution of the galaxy satellite fraction. The evo- 
lution was obtained from data observed in the same rest- 
frame band, and provides for simpler interpretations as com- 
pared to previous studies using data from different restframe 
bands. Various luminosity threshold samples at different red- 
shifts were selected and the corresponding best-fit HOD pa- 
rameters for two similar HOD models obtained. This is done 
in order to single out possible degeneracies and inconsis- 
tencies with the fitting procedure at high redshifts. On the 
whole, both models are in agreement with each other and 
show similar trends in evolution. The impact of our selection 
on the average halo mass is addressed using the Millennium 
simulation. We find that a growth in halo mass as seen in the 
data could rather be an underestimation of ~ 10% to what 
is seen in an 'ideal' sample containing all the descendants. 
Therefore a measure in the growth of mass of a halo can be 
mainly attributed to the hierarchical formation of structure 
and not due to the typology of the selection. 

We find that the number-weighted average halo mass 
grows by ~ 90% from redshift 1.0 to 0.5. This is the first 
time a growth in the underlying halo mass has been mea- 
sured at high redshifts within a single data survey, and 
provides evidence for the rapid accretion phase of massive 
halos. The mass accretion history follows the form given 



in Wechsler et al. (2002) with M{z) = Moe"P^ where 
/3 ~ 1.07 ± 0.57(1.54 ± 0.13) when only the VVDS points 
were used and /3 ~ 1.94 ± 0.10(2.09 ± 0.04) after includ- 
ing the SDSS data as reference points at low redshift and 
depending on the model used to obtain the best fits. The 
addition of the low redshift SDSS points adds complications 
due to the addition of possible systematics by comparing 
data from two different rest-frame bands, even after conver- 
sion to a common fiducial band. We adopt the average value 
of /3 ~ 1.3 ± 0.30 from the VVDS points when discussing a 
growth in halo mass, and found to be slightly higher than 
the results from N-body simulations. 

If we express this result in terms of the expected halo 
mass at present times, Mq ~ 10^^'^ Mq, such halos ap- 
pear to accrete m ~ 0.25 Afo between redshifts of 0.5 and 
1.0. Stewart et al. (2007) have shown that ~ 25% (80%) 
of Mq = lQ^^h~^ Mq halos experienced an m > 0.3A/o 
(m > Q.lMo) merger event in the last 10 Gyr, this would 
translate into a m > O.lAfo merger event over the redshift 
range z= [0.5- 1.0] for the high mass halos here. From merger 
rate studies one finds that 30% of the stellar mass of massive 
galaxies with lO^'^M© < M < lO^M© has been assembled 
through mergers since z=l (de Ravel et al. 2009, and refer- 
ences therein). The integrated stellar mass growth obtained 
can then be compared to the halo mass growth obtained 
here. 

For samples at similar redshifts we see that the average 
halo mass, < M >, generally increases with the luminosity 
threshold of the sample, with a very mild hint of a decreasing 
galaxy satellite fraction. This implies that galaxies in the 
faint sample show a stronger probability of being satellites 
in low mass halos as compared to bright galaxies in massive 
halos. 
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Figure 6. The best fit masses (Mmin, Mi) and power law slope a of halos for Z model versus (Lt^i^e^^/L, ). The different symbols 
represent results for the VVDS, where each set of symbols correspond to samples lying in different redshift ranges (circles, solid triangles, 
and empty squares represent samples having z = [0.4 — 0.6], [0.6 — 0.8], [0.8 — 1.1] respectively). The superimposed dashed lines depict 
the SDSS. 



We also find that the satellite fraction or average num- 
ber of satellite galaxies appears to slowly increase over the 
redshift interval [0.5,1.0], but a stronger increase by a factor 
of ~3 over z=[0. 1,0.5] is seen. This can be understood in 
terms of the dynamical friction that subhalos hosting satel- 
lite galaxies encounter within their host halos. The efficiency 
of dynamical friction depends on the relative subhalo to halo 
mass. Subhalos experience more efficient dynamical friction 
in low mass halos, which can be thought of as progenitor 
halos at high redshift. The subhalos are continuously sub- 
jected to tidal stripping and gravitational heating within 
the dense environments and get eroded if not completely. 
As time evolves the halo accretes mass and undergoes merg- 
ers with other halos. The subhalos that form as remnants 
of halo mergers are now more likely to remain intact within 
the higher mass halo, in turn leading to a larger number of 
satellite galaxies in present-day halos. 



A comparison with the SDSS results shows a few inter- 
esting features. The value for Mmin, which is the mass of a 
halo hosting at least one central galaxy on average, in 40% of 
the luminosity threshold VVDS samples is similar to values 
for local SDSS galaxies. Whereas, M\ is generally higher for 
VVDS galaxies as compared to what is seen locally. The ra- 
tio of M\/Mmin is fouud to be considerably higher (almost 
a factor of 2) in the VVDS as compared to the SDSS results. 
This shows that in order to begin hosting satellite galaxies, 
halos at high redshift need to accrete a larger amount of 
mass than is seen locally. Hence one would observe roughly 
twice as many local satellite galaxies than high redshift ones 
within the same evolved halo mass. This is another line of 
evidence in favor of the lower observed satellite fraction at 
high redshift and high local satellite fraction. This interpre- 
tation is highly simplified in light of the fact that the results 
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have been obtained with data taken in different restframe 
bands. 

In order to investigate further and better constrain the 
mass growth and evolution in the number of satelMte galax- 
ies per halo over a larger rcdshift range, one needs to have 
samples from the same survey at low redshifts. This can be 
done with samples from deeper and wider redshift surveys. 
Here we have concentrated on luminosity-threshold samples 
leading to a link between the luminosity of galaxies and the 
underlying dark matter distribution. The present paper can 
be seen as a precursor to many studies that can be carried 
out with larger samples than the VVDS, including CLF - 
conditional luminosity function studies (e.g. van den Bosch 
et al. 2003, etc.), analyses with galaxy samples of different 
stellar masses (Zheng et al. 2007), etc.. They will certainly 
add to the understanding of the vast pool of underlying dark 
matter properties and hopefully obtain tighter constraints 
on models of galaxy formation. 



6 ACKNOWLEDGEMENTS 

UA would like to acknowledge funding from the Mario 
Curie training network supported by the European Commu- 
nity's Sixth Framework Programme (FP6). UA also thanks 
Alessandro Sozzetti and Martin Kilbinger for helpful dis- 
cussions. The authors thank the anonymous referee for very 
useful comments and suggestions that helped improve the 
paper. This research program has been developed within 
the framework of the VVDS consortium. This work has 
been funded in part by the ANR program ANR-05-BLAN- 
0283 and partially supported by the CNRS-INSU and its 
Programme National de Cosmologie (France), and by the 
Italian Ministry (MIUR) grants COFIN2000 (MM02037133) 
and COFIN2003 (No. 2003020150) and by INAF grants 
(PRIN-INAF 2005) and the grant of Polish Ministry of Sci- 
ence and Higher Education PBZ/MNiSW/07/2006/34. The 
VLT-VIRMOS observations have been carried out on guar- 
anteed time (GTO) allocated by the European Southern 
Observatory (ESO) to the VIRMOS consortium, under a 
contractual agreement between the Centre National dc la 
Recherche Scientifique of France, heading a consortium of 
French and Italian institutes, and ESO, to design, manufac- 
ture and test the VIMOS instrument. 



REFERENCES 

Benson A., Cole S., Frenk C.S., Baugh CM., Lacey C.G., 

2000, MNRAS, 311, 793 
Berhnd A., Weinberg D.H., 2002, ApJ, 575, 587 
Blaizot J., et al., 2005, MNRAS, 360, 159 
Carroll S.M, Press W.H., Turner E.L., 1992, ARA&A, 30, 

499 

Coil A., et al., 2006, ApJ, 644, 671 
Connolly A., et al. 2002, ApJ, 579, 42 
Conroy C, Wechsler R.H., Kravtsov A.V., 2006, ApJ, 647, 
201 

Conroy C, et al. 2007, ApJ, 654, 153 

Cooray A., Sheth R.K., 2002, Phys Reps, 372, 1 

Cooray A., 2006, MNRAS, 365, 842 

Davis M., Peebles P.J.E., 1983, ApJ, 267, 465 



de Ravel L. et al., 2009, A&A, 498, 379 

Gott J.R., Turner E.L., 1979, ApJ, 232, L79 

Gunn J.E., Gott J.R.III, 1972, ApJ, 176, 1 

Guzzo, L., lovino. A, Chincarini, G., Giovanelli, R. & 

Haynes, M.P., 1991, ApJ, 382, L5 
Hamana T., Ouchi M., Shimasaku K., Kayo I., Suto Y., 

2004, MNRAS, 347, 813 
Hawkins E., Maddox S., Cole S., et al. (the 2dFGRS Team) 

2003, MNRAS, 346, 78 
Ilbert O., et al. 2005, A&A, 439, 863 
Jenkins A. et al., 1998, ApJ, 499, 20 

Kauffmann G., Colbcrg J.M., Diaferio A., White S.D.M., 

1999, MNRAS, 303, 188 
Kravtsov A.V., Berlind A. A., Wechsler R.H., Klypin A. A., 

Gottlober S., AUgood B., Primack J.R., 2004, ApJ, 609, 

35 

Landy S.D., Szalay A.S., 1993, ApJ, 412, 64 

Le Fevre O., et al., 2003, The Messenger, 111, 18 

Le Fevre O., et al., 2004, A&A, 428, 1043 

Le Fevre O., et al., 2005, A&A, 439, 845 

Le Fevre O., et al., 2005, A&A, 439, 877 

Lee K.-S., Giavalisco M., Gnedin O.Y., Somerville R.S., 

Ferguson H.C., Dickinson M., Ouchi M., 2006, ApJ, 642, 

63 

Li Y., Mo H.J., van den Bosch E.G., Lin W.P., 2007, MN- 
RAS, 379, 689 

Magliocchetti M., Porciani C, 2003, MNRAS, 346, 186 

Markwardt C., 2009, ASPC, 411, 251 

Mc Crackcn H., et al., 2003, A&A, 410, 17 

Meneux B., et al. 2006, A&A, 452, 387 

Monaco P., Theuns T., Taffoni G, 2002, MNRAS, 331, 587 

More J., Numerical Analysis, ed. G. A. Watson (Springer- 
Verlag: Berlin), 630, 105 

Navarro J., Frenk C, White S.D.M., 1997, ApJ, 490, 493 

Ouchi M., Hamana T., Shimasaku K., et al. 2005, ApJL, 
635, 117 

Peacock J. A., 1997, MNRAS, 284, 885 
Peebles P.J.E., 1974, A&A, 32, 197 

Phleps S., Peacock J.A., Meisenheimer K., Wolf G., 2006, 

A&A, 457, 145 
Polio A. et al., 2005, A& A, 439, 887 
Polio A. et al., 2006, A& A, 451, 409 

Press W.H., Teukolsky S.A., Vetterling W.T., Flannery 
B.P., 1992, Numerical Recipes in C, Cambridge University 

Press 

Sheth R.K., Tormen G., 1999, MNRAS, 308, 119 
Springel V., 2005, Nature, 435, 629 

Stewart K., Bullock J.S., Wechsler R.H., Mailer A.H., Zent- 

ner A.R., 2007, astro-ph/0711.5027 
Tinker J.L., Weinberg D.H., Zheng Z., Zehavi I., 2005, ApJ, 

631, 41 

Totsuji H., Kihara T., 1969, PASJ, 21, 221 
van den Bosch F.C., Yang X., Mo H.J., 2003, MNRAS, 340, 
771 

Wechsler R.H., Bullock J.S., Primack J.R., Kravtsov A.V., 

Dekel A., 2002, ApJ, 568, 52 
Weinberg D.H., Dave R., Katz N., Hernquist L., 2004, ApJ, 

601, 1 

White S.D.M., Rees M.J., 1978, MNRAS, 183, 341 
Yang X., Mo H.J., van den Bosch F.C., 2003, MNRAS, 339, 
1057 



VVDS: Evolution m the HON over z=0-1.3 13 



Zehavi I., Weinberg D.H., Zheng Z., et al. 2004, ApJ, 608, 
16 

Zehavi I., et al., 2005, ApJ, 630, 1 

Zheng Z., 2004, ApJ, 610, 61 

Zheng Z., Coil A., Zehavi I., 2007, ApJ, 667, 760 



List of Affiliations 

^Laboratoire d'Astrophysique de Marseille, UMR 6110 
CNRS-Universite de Provence, BPS, F-13376 Marseille 
Cedex 12, France 

^INAF-Osservatorio Astronomico di Brera, Via Brera 28, 
1-20021, Milan, Italy 

^Max Planck Institut liir Extraterrestrische Physik 
(MPE), Giessenbachstrasse 1, D-85748 Garching bei 
Miinchen , Germany 

*Max Planck Institut liir Astrophysik, D-85741, Garching, 
Germany 

''Centre de Physique Theorique, UMR 6207 CNRS- 
Universite de Provence, F-13288, Marseille, France 
^Universitatssternwarte Miinchen, Scheinerstrasse 1, D- 
81679 Miinchen, Germany 

^The Andrzej Soltan Institute for Nuclear Studies, ul. Hoza 
69, 00-681 Warszawa, Poland 

^Astronomical Observatory of the Jagiellonian University, 
ul Orla 171, PL-30-244, Krakow, Poland 
^INAF-Osservatorio Astronomico di Bologna, Via Ranzani 
1, 1-40127, Bologna, Italy 

i°IASF-INAF, Via Bassini 15, 1-20133, Milano, Italy 
"IRA-INAF, Via Gobetti 101, 1-40129, Bologna, Italy 
^■^INAF-Osservatorio Astronomico di Capodimonte, Via 
Moiariello 16, 1-80131, NapoU, Italy 

^^Universita di Bologna, Dipartimento di Astronomia, Via 
Ranzani 1, 1-40127, Bologna, Italy 

^'*Laboratoire d'Astrophysique de Toulouse/Tabres 
(UMR5572), CNRS, Universite Paul Sabatier - Toulouse 
III, 

Observatoire Midi-Pyrenees, 14 av. E. Belin, F-31400, 
Toulouse, France 

^''Institut d'Astrophysique de Paris, UMR 7095, 98 bis Bvd 
Arago, F-75014, Paris, France 

^^Observatoire de Paris, LERMA, 61 Avenue de 
rObservatoire, F-75014, Paris, France 

^^Astrophysical Institute Potsdam, An der Sternwarte 16, 
D-14482, Potsdam, Germany 

^^INAF-Osservatorio Astronomico di Roma, Via di Frascati 

33, 1-00040, Monte Porzio Catone, Italy 

^^Universita di Milano-Bicocca, Dipartimento di Fisica, 

Piazza delle Scienze 3, 1-20126, Milano, Italy 

^"Integral Science Data Centre, ch. d'Ecogia 16, CH-1290, 

Versoix, Switzerland 

^^Geneva Observatory, ch. des Maillettes 51, CH-1290, 
Sauverny, Switzerland 

^■^Centro de Astroffsica da Universidade do Porto, Rua das 
Estrelas, P-4150-762, Porto, Portugal 

■^•^Institute for Astronomy, 2680 Woodlawn Dr., University 
of Hawaii, Honolulu, Hawaii, 96822, USA 
^''School of Physics & Astronomy, University of Nottingham 
University Park, Nottingham, NG72RD, UK 
Canada France Hawaii Telescope corporation, Mamalahoa 



Hwy, Kamuela, Hl-96743, USA 

^^INAF - Osservatorio Astronomico di Torino, Strada 
Osservatorio 20, Pino Torinese, 1-10025, Italy 



7 APPENDIX 

Here we present the table of values for the projected cor- 
relation function (wp) and associated errors (cr) at different 
values of Vp (in units of Mpc hT^) for the different subsam- 
ples mentioned in Table[T]in the main text. The values for Wp 
and a are reported horizontally at the corresponding values 
for r-p in the top row. 
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Table 5. Projected correlation function with associated errors at different rp for the different subsamples 
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